Nonclassical effects in cold trapped ions inside a cavity 
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We investigate the dynamics of a cold trapped ion coupled to the quantized field inside a high- 
finesse cavity, considering exact resonance between the ionic internal levels and the field (carrier 
transition) . We derive an intensity-dependent hamiltonian in which terms proportional to the square 
of the Lamb-Dicke parameter (rj) are retained. We show that different nonclassical effects arise in 
the dynamics of the ionic population inversion, depending on the initial states of the vibrational 
motion/field and on the values of rj. 
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I. INTRODUCTION 



The manipulation of simple quantum systems such as trapped ions pi has opened new possibilities regarding not 
only the investigation of foundations of quantum mechanics, but also applications on quantum information. In such 
a system, internal degrees of freedom of an atomic ion may be coupled to the electromagnetic field as well as to the 
motional degrees of freedom of the ion's center of mass. Under certain circumstances, in which full quantization of the 
three sub-systems becomes necessary, we have an interesting combination of two bosonic systems coupled to a spin 
| like system. A possibility is to place the ion inside a high finesse cavity in such a way that the quantized field gets 
coupled to the atom. A few papers discussing that arrangement may be found in the literature, e.g., the investigation 
of the influence of the field statistics on the ion dynamics |^|, |) , the transfer of coherence between the motional states 
0^ ! and the field , a scheme for generation of matter-field Bell-type states || , and even propositions of quantum logic 
gates ||. Several interaction hamiltonians analogue to the ones found in quantum optical resonance models may be 
constructed even in the case where the field is not quantized, but having the center-of-mass vibrational motion playing 
the role of the field. For instance, if we consider a two-level atom having atomic transition too and center-of-mass 
oscillation frequency of v in interaction with a laser of frequency lu^, if lo^ — ujq = ~ v (laser tuned to the first red 
Oh, sideband), it results an interaction hamiltonian of the Jaynes-Cummings type 

c"| ■ Hi = ih£l(<7+a — o-d)). (1) 

If lol — loq = v (laser tuned to the first blue sideband), the resulting hamiltonian is an "anti- Jaynes-Cummings" type, 

> '■ H l = ihn((T + a< - cr-a). (2) 

• <-H , 

Note that here the bosonic operators, a^,a, are relative to the center-of-mass oscillation motion. The ion itself is 
considered to be confined in a region much smaller than the laser light wavelength (Lamb-Dicke regime). Other forms 
of interaction hamiltonians may be constructed in a similar way U ||] . 

In this paper we explore further the consequences of having the trapped ion in interaction with a quantized field. 
We show that under certain conditions, the quantum nature of the field is able to induce intensity-dependent effects in 
the trapped ion dynamics, somehow analogue to models in cavity quantum electrodynamics . We shall remark 

that here we derive a hamiltonian with an intensity-dependent coupling from a more general hamiltonian, which is 
different from the phenomenological approach discussed in ^] . 



II. THE MODEL 



We consider a single trapped ion, within a Paul trap, placed inside a high finesse cavity, and having a cavity mode 
coupled to the atomic ion. The vibrational motion is also coupled to the field as well as to the ionic internal degrees 
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of freedom, in such a way that the hamiltonian will read |3[ ^). 



H = hva)a + huii^b + h^-a z + hg(a + + (T-)(P + b) cos r;(a^ + a) 



(3) 



where (a) denote the creation (annihilation) operators of the center-of-mass vibrational motion of the ion (frequency 
v) , w (b) are the creation (annihilation) operators of photons in the field mode (frequency us) , uiq is the atomic frequency 
transition, g is the ion-field coupling constant, and rj = 27rao/A is the Lamb-Dicke parameter, being ao the amplitude 
of the harmonic motion and A the wavelength of light. Note that the ion's position inside the cavity is different 
from the one considered in 0]. Tipically the ion is well localized, confined in a region much smaller than light's 
wavelenght, or rj <c 1 (Lamb-Dicke regime). Usually expansions up to the first order in rj are made in order simplify 
hamiltonians involving trapped ions, which results in Jaynes-Cummings like hamiltonians such as the ones in Eqs. d3) 
and (|^). However, even for small values of the Lamb-Dicke parameter, we show that in the situation here discussed, 
an expansion up to second order in r/ becomes necessary. Several interesting effects, such as long time scale revivals 
will depend on those terms, as we are going to show. We may expand the cosine in Eq. (|^) and obtain 



cos rj(a) + a) « 1 



77 2 (l + 2ata) rf{a) 2 + a 2 ) 



(4) 



The interaction hamiltonian will become 

Hi = fig 



7? 2 (l + 2qta) ^(a^ + a 2 ) 
2 2 



(a + +a_)(6 t +6). 



(5) 



We may then rewrite the hamiltonian above in the interaction picture, in order to apply a rotating wave approximation. 
If we tune the light field so that it exactly matches the atomic transition, i.e., ujq — lu = (carrier transition), and 
after discarding the rapidly oscillating terms, we obtain the following (interaction picture) interaction hamiltonian: 



H! = hg 



1 



77 2 (l + 2a t fi) 



(CT_6 f +a + b). 



(6) 



The resulting hamiltonian is alike a Jaynes-Cummings hamiltonian. It describes the annihilation of a photon and a 
simultaneous atomic excitation (w and b are field operators) , but having an effective coupling constant which depends 
on the excitation number (era) of the ionic oscillator. As we are going to show, this will bring interesting consequences 
on the ion's center-of-mass dynamics. We would like to point out that we have retained terms of the order of rj 2 in 
the co-sine expansion, which are much smaller than one. However, the product r] 2 (a^a} will not be negligible for a 
sufficiently large excitation number of the vibrational motion. 

The evolution operator associated to the hamiltonian (^) is given by 



U{t) = C m - n+ i\e)(e\ + C m - n \g)(g\ - iS m - n+1 b\e)(g\ - itfS., 



m;n+l 



\g)< 



where 



a 



m;n+l 



COS 



(.g[l-r ? 2 (l + 2m)/2] y/(n + l)t 



(7) 



(8) 



Cm;n = cos (g [l - rj 2 (1 + 2m)/2] Vfit) , (9) 



and 



sin 



' K g [1 - 7y 2 (l + 2m)/2] y/(n + l) tj 

S m -n+l = 7= , (10) 

V(™ + 1) 



where we have used the notation rh = a) a and h = wb. 
We consider now the following (product) initial state 



P(0) = |e><e| <8> p/(0) <8) p„(0), 



(11) 
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or ion's internal levels prepared in the excited state |e)(e|, the field prepared in a generic state /5/(0), and the ion's 
vibrational center of mass motion prepared in a state p v (0). Its time evolution, governed by the evolution operator 
(0) results, at a time t, in the following (joint) state, 

p(t) = <?m;n+iP/(0)/8 u (0)C min+ i|e)(e| + b^S m . n+ ipf(0)p v {0)S m . n+1 b\g}(g\ 

+ iC m;n +iPf(0)pv(0)Sm;n+ib\e)(g\ - ^S m . n+1 pf(0)p v {0)C m . n+1 \g)(e\, (12) 

where the operators C and S are the ones in Eq. (|^), (^) and ([lO]) above. 

The state in Eq. ([l2|) is, in general, an entangled state involving the ion's internal (electronic) degrees of freedom, 
the vibrational motion as well as the cavity field. We will now investigate the atomic population inversion considering 
simple initial states for the center-of-mass vibrational motion and the cavity field. 



III. ATOMIC DYNAMICS 



The (internal level) ionic dynamics will depend on the distributions of initial excitations of both the field and the 
center center-of-mass vibrational motion, given by (n|p/(0)|n) = p^.„(0), (m\p v (Q)\m) — Pm m (0) respectively. For 
instance, the atomic population inversion may be written as 

oo 

W(t) = Tr [a z p(t)} = E pLMp"m,m(0) cos [2g (l - if {I + 2m)/2) vVTTt] . (13) 

n,m— 

Due to the different frequencies (y — 2g[l — i] 2 (l + 2m)]y/n + 1) in the co-sine argument, we expect different 
structures of beats in the inversion due to different preparations of the center-of-mass motion (index m, without a 
square root) and the field (index n, with y/n+1). We immediately identify two well known particular cases - if the 
vibrational motion is prepared in a number state (Pm. m (0) = Sk. m ), the atomic inversion in Eq. ( |l3| ) will reduce 
to the characteristic pattern of the Jaynes-Cummings model - if otherwise the field is prepared in a number state 
(Pn n(0) = fil,n)) we have a periodic behavior for the atomic inversion (no square root in the co-sine argument). This 
is shown in Fig. 1, where the field is initially in the vacuum state and the ion in a coherent state of the center-of-mass 
motion. The Lamb-Dicke parameter r\ is typically less than unity, and the lowest order term in the expansion above, 
[see Eq.(^)] is of second order in r\. Nevertheless, as we are going to show, interesting effects will arise if the term 
r/ 2 (l + 2m) becomes large enough Q 

We now consider both the center-of-mass motion and the field initially prepared in coherent states and the atom 
in the excited state |e). After performing the summation over m in Eq. (|l3]) , we obtain the following expression 

00 

W(t) = Y,Pn,n(0)Wn(t), (14) 

with 



w n (t) = exp {— to [l — cos (277 2 v / n + 1 gt)] } 

x cos |2 ^1 - y^j V" + 1 gt - m sin [2r] 2 Vn + 1 gt] X . (15) 



The dynamics of the population inversion predicted by Eq. (|13|) may be interpreted in terms of two families of 
"revival" times. The revival times associated to the field (t(), when terms n and n + 1 arc in phase Jll| 

f 2k^\ffi , 
g*r ~ „2 (1+2 -=y ; fc=l,2,..., (16) 

1 2 

depend on m, the mean excitation number of the center-of-mass motion, and on n, the mean excitation number of 
the field. On the other hand, the revival times associated to the vibrational motion (£"), when terms m and m + 1 
are in phase, depend only on n, the mean excitation number of the field, or 

gtl « k Z- J ; k=l,2,..., (17) 
r] z \>n + 1 
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It is also possible to estimate collapse times in either case, or 




1 



(18) 



. 2 (l+2m) ' 
2 



and 



i 



(19) 



j} 2 m(n + 1) 



We note that distinct patterns for the atomic inversion will arise depending on the values of rj 2 , m and n, which are 
quantities determining the revival/collapse times. In order to appreciate those different situations, we show, in what 
follows, some plots of the atomic inversion for different excitations of the ionic motion and held as well as for different 
values of 77. 

In Fig. 2 we have a plot of the atomic inversion as a function of time, having r/ — 0.02, m = 4 and n = 25. It is 
shown the usual pattern of collapses and revivals. In this case the relevant times are so that gt{ w 1 and gtl w 31. 
For a larger Lamb-Dicke parameter, 77 = 0.04, the oscillations in the atomic inversion are attenuated, as we see in 
Fig. 3. This means that one has to be careful while expanding over the parameter 77, since in the model we are here 
discussing, terms of the order of rj 2 have a significant influence on the system evolution. Another interesting behavior 
predicted by Eq. ( fL3| ) is the phenomenon of "super-revivals" of the atomic inversion, or revivals taking place at long 
times. Those may occur in driven atoms | fLl) or in trapped ions under certain conditions Jl^] . Our model also predicts 
"super revivals" if both the center-of-mass motion and the field are prepared in coherent states, as it is shown in Fig. 
4 (same parameters as in Fig. 2). The "super revival" time in this case is given by gt v r w 1540, with a collapse time 
gt v c ~ 245. Still for long times, if 77 = 0.04, as it is shown in Fig. 5, there will be changes in the "super revival" time, 
gt v r « 385 as well as in the collapse time, gt v c ~ 61. In this case the revival time for "ordinary revivals" is very close to 
the collapse time associated to the long time scale dynamics, so that the revivals themselves are attenuated (compare 
Fig. 2 with Fig. 3). Different patterns will occur for different values of m and n. For instance, if m = 25, n — 4 and 
rj = 0.2, we have gt v r « 48 and gtl sa 14. In this case there is only one oscillation in the collapse region, as shown in 
Fig. 6. A different pattern for the oscillations occurs if the collapse time gt v c ~ 24 is slightly less than the revival time 
gt' r w 26 (see Fig. 7). The short time revivals are strongly suppressed, as we see in Fig. 7. We may compare the 
atomic inversion in Fig. 7 to the one in Fig. 2; in the latter case the collapse time is long enough to still allow revivals, 
while in the former one a shorter collapse time significantly reduces the amplitude of the revivals (at gt( «j 26). Apart 
from that, we also note a modulation in the oscillations, which is basically due to the combination of oscillations 
originated from the field and vibrational motion quanta distributions. As a final comment, we have that in the case 
of the periodic dynamics (Fig. 1), gtl — and gt v r « 1257, which is in agreement with the revival times shown in Fig. 
1. We would like also to point out that dissipation in the cavity and loss of coherence in the vibrational motion will 
certainly have a destructive action, and effects occurring at longer time scales may not be apparent. However, the 
atomic inversion is highly sensitive to the initial state preparation and the Lamb-Dicke parameter also for short times, 
and the observation of at least some of the effects above described might become possible with further improvements 
in the quality of cavities as well as in the control of trapped ions |13j . 



We have investigated the dynamics of a single trapped ion enclosed in a cavity. We have found that, in the case 
of exact atom-field resonance (carrier transition), terms of the order of rj 2 can make a significant contribution and 
should therefore be retained. This gives rise to a hamiltonian containing an intensity-dependent coupling. Besides, 
the quantum nature of the field plays an important role. We have shown that, the atomic inversion as a function of 
time may display different structures of beats depending on the initial preparation of the electromagnetic field and 
the ionic motion as well as on the Lamb-Dicke parameter 77, and effects such as suppression/attenuation of the Rabi 
oscillations, long time scale revivals as well as a periodic dynamics may occur. Those interesting features may be 
understood in terms of the two collapse times and two revival times characteristic of the dynamics. 



IV. CONCLUSION 
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FIG. 1: Atomic inversion as a function of the scaled time gt. The field is initially in the vacuum state, n = 0, and the vibrational 
motion in a coherent state with m = 4. The ion is prepared in its internal excited state, and rj = 0.05. 
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FIG. 2: Atomic inversion as a function of the scaled time gt. The field is initially in a coherent state, n = 25, and the vibrational 
motion in a coherent state with m = 4. The ion is prepared in its internal excited state, and r\ = 0.02. 
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FIG. 3: Atomic inversion as a function of the scaled time gt. The field is initially in a coherent state, n = 25, and the vibrational 
motion in a coherent state with m = 4. The ion is prepared in its internal excited state, and rj = 0.04. 
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FIG. 4: Atomic inversion as a function of the scaled time gt. The field is initially in a coherent state, n = 25, and the vibrational 
motion in a coherent state with m = 4. The ion is prepared in its internal excited state, and r\ = 0.02. 
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FIG. 5: Atomic inversion as a function of the scaled time gt. The field is initially in a coherent state, n = 25, and the vibrational 
motion in a coherent state with rn — 4. The ion is prepared in its internal excited state, and rj = 0.04. 
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FIG. 6: Atomic inversion as a function of the scaled time gt. The field is initially in a coherent state, n = 1.69, and the 
vibrational motion in a coherent state with m = 10.24. The ion is prepared in its internal excited state, and rj = 0.2. 
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FIG. 7: Atomic inversion as a function of the scaled time gt. The field is initially in a coherent state, n = 16, and the vibrational 
motion in a coherent state with m = 16. The ion is prepared in its internal excited state, and r\ = 0.05. 



